###########################################
# Continuos exposure 
###########################################

rm (list = ls())

library(Hmisc)
library(readstata13)
library(foreign)
library(lmtest)
library(sandwich)
library(stargazer)
library(msm)
library(ggplot2)
library(did)
library(lfe)
library(gridExtra)

# load data
load("one_week_protests_riots_2024March1.RData")

##############
# Table A8
##############

m1 = felm(support_eviction  ~ t_con + year_f + t_con*year_f + woman + age + education | 0 | 0 | idencuesta, data = d, keepCX = TRUE)
t1 = m1
m1 = tidy(m1,conf.int = T)
m1

##############
# Table A9
##############

m2 = felm(support_students  ~ t_con + year_f + t_con*year_f + woman + age + education | 0 | 0 | idencuesta, data = d, keepCX = TRUE)
t2 = m2
m2 = tidy(m2,conf.int = T)
m2
